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Magnetic and chiral universality classes in a 3D Yukawa model* 
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The 3D Yukawa model with U(l) chiral symmetry is investigated in a broad interval of parameters using the 
Binder method. Critical exponents of the Wilson-Fisher (magnetic) and Gross-Neveu (chiral) universality classes 
are measured. The model is dominated by the chiral universality class. However at weak coupling we observe a 
crossover between both classes, manifested by difficulties with the Binder method which otherwise works well. 



1. Introduction 

The existence of nontrivial fixed points in 4D 
is not yet ruled out in nonperturbative calcula- 
tions. Several gauge theories have, in addition 
to the Gaussian fixed points, also suspicious crit- 
ical points at strong coupling. One of the pur- 
poses of our investigation of the 3D Yukawa model 
(Y3) was to learn how to deal with models hav- 
ing several nontrivial fixed points and crossovers 
between the corresponding universality classes. 

The Y3 model is known to have two nontriv- 
ial fixed points: the Wilson-Fisher fixed point of 
the pure scalar 0| theory at vanishing Yukawa 
coupling with a magnetic type phase transition, 
and the fixed point of the 3D Gross-Neveu model 
(GN3) with a chiral phase transition. 

The model is also interesting from the point 
of view of statistical mechanics. Transitions be- 
tween different universality classes have been in- 
vestigated in spin models but not yet in 
fermionic ones. The most promising method used 
is Binder's method of finite size scaling analysis. 
We have applied this method succesfuUy to the 
chiral phase transition. The failure of this method 
also indicates the occurence of a crossover to the 
magnetic universality class. The existence of in- 
termediate universality classes between the two, 
corresponding to the known fixed points, is also of 
interest [|l|, but in the Y3 model we didn't detect 
any signs for it. 



*Work supported by the DFG and BMBF. Computations 
have been performed on the CRAY-YMP in Jiilich and on 
the Quadrics QH2 in Bielefeld. 



2. Lattice action and phase diagram 

We studied the Y3 model on the lattice with 
staggered fermions, hypercubic Yukawa coupling 
and U(l) chiral symmetry. The action is 



S = 



E 



-2fi 



2 E/ ^^'f^ 



(1) 



The staggered sign factors are rj^ i — 1, ij^.f^ ~ 
(-l)^i+-+^^'-i and £^ = (-1)^1+'-+^^. The in- 
dices X + fi and X + b denote, respectively, the 
nearest neighbors of the site x and the corners of 
the associated elementary cube. The scalar field 
is a two component real field. 
From considerations of the effective potential to 
1-loop order and numerical simulations we deter- 
mined the phase diagram sketched in fig. |l|. It has 
two well-known limits: the (/)| model a.t y — and 
the GN3 model at A = 0, k = 0. The region be- 
low the critical surface is the paramagnetic phase 
(PM), the region above the ferromagnetic phase 
(FM). For negative values of the parameter k we 
further expect an antiferromagnetic phase (AFM) 
and a ferrimagnetic phase (FI). 

3. Numerical methods 

For the numerical simulation we used a hy- 
brid Monte-Carlo program. Ferrenberg-Swendsen 
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Figure 1. Sketch of the Y3 phase diagram. 



multihistogram reweighting was used to interpo- 
late between the measured points in the parame- 
ter space. 

In order to find an appropriate method to iden- 
tify the universahty class corresponding to a cer- 
tain critical point we tried several methods to 
compute critical exponents in the pure scalar 
limit of the model (y = 0). The direct method us- 
ing the scaling laws of the physical quantities and 
the method using the Lee- Yang zeros of the parti- 
tion function needed significantly more statistics 
compared to Binder's method of finite size scal- 
ing analysis of cumulants. Thus we used only the 
Binder method in the fermionic case, with the 
following definition for the cumulant: 



Xi,X2 i 



(2) 



The validity of the hyperscaling hypothesis im- 
plies that these cumulants are independent of the 
lattice size at the critical point. This is a very 
precise way of determining critical couplings. 

Critical exponents can be computed by consid- 
ering pairs (6L, L) of lattice sizes. The cumulants 
deliver the exponent v of the scalar correlation 
length ^: 

dUbL 
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dUr 



V 



(3) 



Similarly we obtained ^ /v and fi/v from the sus- 
ceptibility XL and the magnetization M^. 



4. The 



b\ model 



The action (|l|) describes at vanishing y free 
massless fermions and 0(2) invariant scalars with 
selfinteraction (0| model). We have determined 
the phase diagram of this model (for the data 
see H]). At positive n a critical line of second 
order phase transitions Kc(A) separates the PM 
and FM phases. At the point Kc(A = 0) the 
theory is dominated by a Gaussian fixed point, 
being asymptotically free at large momenta. All 
other points on the critical line /«c(A > 0) lead to 
a nontrivial continuum limit, dominated by the 
Wilson-Fisher fixed point. In order to test this 
expectation, we measured the renormalized quar- 
tic coupling and the critical exponents v, (3/v and 
at two values of the bare scalar quartic cou- 
pling: A = 0.5 and A = 00. Fig. 3a illustrates 
the quality of determination of Kc by the Binder 
method at A = 00. These measurements helped 
us to develop and test the methods later used for 
nonvanishing y. 

The renormalized scalar quartic coupling 



A_R = {LamjifUL 



(4) 



has been determined by keeping the ratio L/^ 
fixed to 4 and varying the lattice size L between 
6 and 12. The extrapolation to infinite volume 
suggests Afl — 26(4) for both values of A. 

The critical exponents we calculated by the 
Binder method are summarized in the following 
table. 



A 




V 


(ijv 




00 
0.5 


0.2275(10) 
0.241(1) 


0.673(19) 
0.687(19) 


0.51(3) 
0.56(5) 


2.03(6) 
1.91(6) 



The results are consistent within error bars and 
also with the hyperscaling hypothesis. They con- 
firm that the continuum limits at A = 00 and 
A = 0.5 belong to the Wilson-Fisher universality 
class. 
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5. Ya model at A = 

The GN3 model with U{1) chiral symmetry 
arises in the hmit A = 0. It is renormahzable 
in the 1/Np expansion and its /3 function has 
been calculated up to 0(l/iV|) |||. The criti- 
cal exponent v resulting from these calculations 
is 1/ ~ 1.0(1). 

The phase structure of the GN3 model can be 
read of fig. |l|. The symmetric phase {y < yc), 
where fermions are massless, is dominated by the 
Gaussian fixed point at y = 0. The critical point 
yc = 1.091(5) is an UV stable nontrivial fixed 
point. 




Figure 2. Determination of the critical exponent 
v in the Gross-Neveu. 



We applied the Binder method to the GN3 
model. We approached the critical point both by 
varying y at A = k = (GN case) and k at A = 
and y — yc- Fig. |^ shows the determination of 1/ 
using eq. (||). Fig. 3b shows that the cumulants 
cross in the GN case at y — 1.091. The obtained 
exponents are perfectly consistent with each other 
and the theoretical values. Their averages are: 
1/ = 1.03(10), P/iy = 0.89(7), j/v = 1.17(9), 
values which are significantly different from the 
Wilson-Fisher exponents and allow the investiga- 
tion of crossover effects between these universality 
classes. 




095 1.100 




Figure 3. Determination of Kc- (a) in the XY3 
model, (b) in the GN3 model and (c) at A = 00, 
y = 0.6. It works well in the first two cases, but 
only poorly in the third case. 

6. Y3 model at A = cxd 

There are theoretical arguments jj] which pre- 
dict that Y3 and GN3 models belong to the 
same universality class for small values of A. We 
tested this conjecture for A = 00 by determining 
the critical exponents at large Yukawa coupling 

(y = i-i)- 

The Binder method works as well as in the GN3 
case. At y = 1.1 and A = 00 we varied k and 
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determined first its critical value. The intersec- 
tion point of the cumulants on different lattice 
sizes delivers with good accuracy Kc — 0.007(2). 
Using eq. (||) the critical exponent 1/ has been 
extracted: i' = 0.89(6). Further we determined: 
P/iy = 0.80(8) and -y/v = 1.30(7). 

All these exponents are consistent with those 
of the GN3 fixed point and significantly different 
from the Wilson-Fisher exponents. We conclude 
that, at least for strong enough Yukawa coupling, 
the critical surface from A = to A = 00 belongs 
to the GN universality class. 

At A = 00 but smaller Yukawa coupling 
crossover phenomena make the determination of 
critical exponents more difficult. 

At y = 0.6 it is difficult to determine Kc (fig- 
3c). For small lattices {L = 6,8,10) the cumu- 
lants cross in the interval k = 0.146 — 0.1466. 
The finite size analysis delivers = 0.79(9), in- 
consistently to the GN universality class. When 
only lattices larger than L = 10 are considered, 
the crossing point of the cumulants is ~ 0.144. 
An analysis for this Kc leads to ly — 0.99(23). 
This value is consistent with the GN critical ex- 
ponent. The reason for the strong dependence of 
the calculated v on Kc is a larger curvature in the 
functions UbL{UL)- Thus the value of the deriva- 
tive dUbb/dUL depends stronger on the value of 
Kc than in the (/)| and strong y cases. 

The computations aX y — 0.3 revealed even 
stronger crossover effects. Though lattices up to 
L — 32 have been used, the value of Kc couldn't be 
determined. A very strong curvature in UbL{UL) 
makes an accurate calculation of impossible on 
such small lattices. 

7. Boson mass in the Y3 model 

In the Gross-Neveu model the bosonic particles 
can be interpreted as fermion-antifermion bound 
states, because the scalar field introduced in the 
action is auxilliary. In the full Yukawa model and 
its purely scalar limit (/)| the (f> field is dynamical. 
We measured the propagator of the two compo- 
nents of the (f) field in momentum space. While in 
the model this propagator can be fitted with 
a free boson Ansatz, its form is very complex at 
nonvanishing Yukawa coupling. We were able to 



fit it with an Ansatz using renormalized lattice 
perturbation theory. 

Fig. ^ shows such a bosonic propagator in the 
Gross-Neveu case and the fit. The characteristic 
form, which is reproduced well by the fit, is the 
same for all A values, including A = 00. Also 
the behavior of the boson mass across the phase 
transition is very similar both at A = and A = 
00. This is another piece of evidence that the 
Y3 and GN3 models are equivalent (belong to the 
same universality class). 




Chl-2 = 0.60(31) 
Moss = 0.384(12) 
Z = 0.8814(93) 
y^R..2= 0.231 16(82) 



2.sum(1-co={p*i)) -> 

Figure 4. Bosonic propagator in the GN3 model 
a,t X — K — 0, y — 1.1. 
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